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\0 ' 1 Introduction 

There are basically two main approaches to define the mark to market of a contingent claim: one relying on 
' the no-arbitrage assumption and the other related to a hedging portfolio, those two approaches converging in 
fT^ ' the specific case of complete markets. In this paper we focus on the hedging approach. A simple introduction 
to the different hedging and pricing models in incomplete markets can be found in chapter 10 of {13j. 
When the market is not complete, it is not possible, in general, to hedge perfectly an option. One has to 
, specify risk criteria, and consider the hedging strategy that minimizes the distance (in terms of the given 
, criteria) between the payoff of the option and the terminal value of the hedging portfolio. In practice the 
price of the option is related to two components: first, the initial-capital value and second the quantitative 
evaluation of the residual risk induced by this imperfect hedging strategy (due to incompleteness). 
Several criteria can be adopted. The aim of super-hedging is to hedge all cases. This approach yields in 
general prices that are too expensive to be realistic [18]. Quantile hedging modifies this approach allowing 
for a limited probability of loss [20]. Indifference utility pricing introduced in [23] defines the price of an 
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option to sell (resp. to buy) as the minimum initial value s.t. the hedging portfolio with the option sold 
(resp. bought) is equivalent (in term of utility) to the initial portfolio. Global quadratic hedging approach 
was developed by M. Schweizer ([35], [10]): the distance defined by the expectation of the square of the 
difference between the hedging portfolio and the payoff is minimized. Then, contrarily to the case of utility 
maximization, in general that approach provides linear prices and hedge ratios with respect to the payoff. 

In this paper, we follow this last approach either to derive the hedging strategy minimizing the global 
quadratic hedging error for a given initial capital, or to derive both the initial capital and the hedging 
strategy minimizing the same error. Both actions are referred to the objective measure. Moreover we also 
derive explicit formulae for the global quadratic hedging error which together with the initial capital allows 
the practitioner to define his option price. 

We spend now some words related to the global quadratic hedging approach which is also called mean- 
variance hedging or global risk minimization. Given a square integrable r.v. H, we say that the pair (Vq, (p) 
is optimal if (c, v) = [Vq, (p) minimizes the functional E — c — vdS^ . The quantity Vq and process (p 
represent the initial capital and the optimal hedging strategy of the contingent claim H. 

Technically speaking, the global risk minimization problem is based on the local risk minimization one 
which is strictly related to the so-called Follmer-Schweizer decomposition (or FS decomposition) of a square 
integrable random variable (representing the contingent claim) with respect to an (J^t)-semimartingale S = 
A/ + A modeling the asset price: M is an (J-"t)-local martingale and ^ is a bounded variation process 
with Aq = 0. Mathematically, the FS decomposition, constitutes the generalization of the martingale 
representation theorem (Kunita-Watanabe representation), which is valid when 5 is a Brownian motion 
or a martingale. Given a square integrable random variable H, the problem consists in expressing H as 
Hq + ^dS + Lt where ^ is predictable and Lt is the terminal value of an orthogonal martingale L to M, 
i.e. the martingale part of S. In the seminal paper [2l], the problem is treated for an underlying process 
S with continuous paths. In the general case, S is said to satisfy the structure condition (SC) if there 
is a predictable process a such that At = J^^ asd{M)s and a^d{M)s < oo a.s. In the sequel, most of 
the contributions were produced in the multidimensional case. Here, for simplicity, we will formulate all the 
results in the one-dimensional case. 

Hq constitutes in fact the initial capital and it is given by the expectation of H under the so called 
variance optimal signed measure (VOM). Hence, in full generality, the initial capital Vq is not guaranteed 
to be an arbitrage-free price. For continuous processes, the variance optimal measure is proved to be non- 
negative under a mild no-arbitrage condition [41]. Aral (|4] and |3]) provides sufficient conditions for the 
variance-optimal martingale measure to be a probability measure, even for discontinuous semimartingales. 
In the framework of FS decomposition, a process which plays a significant role is the so-called mean variance 
trade-off (MVT) process K. This notion is inspired by the theory in discrete time started by [36]; under 
condition (SC), in the continuous time case K is defined as Kt ~ /p a'j.d{M)s, t e [0,r]. In fact, in [35] 
also appear a slight more general condition, called (ESC), together with a corresponding EMVT process; 
we will nevertheless not discuss here further details. If the MVT process is deterministic, [38] solves the 
mean-variance hedging problem and also provides an efficient relation between the solution of the global 
risk minimization problem and the FS decomposition, see Theorem 14. II We remark that, in the continuous 
case, treated by [21], no need of any condition on K is required. It also shows that, for obtaining the 
mentioned relation, previous condition is not far from being optimal. The next important step was done 
in |30| where, under the only condition that K is uniformly bounded, the FS decomposition of any square 
integrable random variable exists, it is unique and the global minimization problem admits a solution. 



2 



More recently has appeared an incredible amount of papers in the framework of global (resp. local) risk 
minimization, so that it is impossible to list all of them and it is beyond our scope. Four significant papers 
containing a good list of references are [12]) [Z], [H] and |43j . 

In this paper, we are not interested in generalizing the conditions under which the FS decomposition 
exists. The present article aims, in the spirit of a simplified Clark-Ocone formula, at providing an explicit 
form for the FS decomposition for a large class of European payoffs H, when the process S is an exponential 
of additive process which is not necessarily a martingale. From a practical point of view, this serves to 
compute efficiently the variance optimal hedging strategy which is directly related to the FS decomposition, 
since the mean-variance trade-off is for that type of processes deterministic. One major idea proposed 
by Hubalek, Kallsen and Krawczyk in [24], in the case where the log price is a Levy process, consists in 
determining an explicit expression for the variance optimal hedging strategy for exponential payoffs and 
then deriving, by linear combination the corresponding optimal strategy for a large class of payoff functions 
(through Laplace type transform) . Using the same idea, this paper extends results of |24| considering prices 
that are exponential of additive processes and contingent claims that are Laplace-Fourier transform of a 
finite measure. In this generalized framework, we could formulate assumptions as general as possible. In 
particular, our results do not require any assumption on the absolute continuity of the cumulant generating 
function of log(S't), thanks to the use of a natural reference variance measure instead of the usual Lebesgue 
measure, see Section In the context of non stationary processes, the idea to represent payoffs functions as 
Laplace transforms was applied by |26| (that we discovered after finishing our paper) to derive explicit pricing 
formulae and by |19| to investigate time inhomogeneous affine processes. However, the |26| generalization 
was limited to additive processes with absolutely continuous characteristics and to the pricing application: 
hedging strategies were not addressed. 

One practical motivation for considering processes with independent and possibly non stationary incre- 
ments came from hedging problems in the electricity market. Because of non-storability of electricity, the 
hedging instrument is in that case, a forward contract with value = e^^'^'^'^'^'^^F^'' — F^"^) where F'^'^ 
is the forward price given at time t < Td for delivery of IMWh at time T^. Hence, the dynamics of the 
underlying 5'° is directly related to the dynamics of forward prices. Now, forward prices are known to exhibit 
both heavy tails (especially on the short term) and a volatility term structure according to the Samuelson 
hypothesis |34| . More precisely, as the delivery date Td approaches, the forward price is more sensitive to 
the information arrival concerning the electricity supply-demand balance for the given delivery date. This 
phenomenon causes great variations in the forward prices close to delivery and then increases the volatility. 
Hence, those features require the use of forward prices models with both non Gaussian and non stationary 
increments in the stream of the model proposed by Benth and Saltyte-Benth, see [3] and also [5]. 

The paper is organized as follows. After this introduction we introduce the notion of FS decomposition 
and describe global risk minimization. Then, we examine at Section 3 the explicit FS decomposition for 
exponential of additive processes. Section 4 is devoted to the solution to the global minimization problem. 
Section 5 to theoretical examples and Section 6 to the case of a model intervening in the electricity market. 
Section 7 is devoted to simulations. 
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2 Preliminaries on additive processes and Follmer-Schweizer de- 
composition 

In the whole paper, T > 0, will be a fixed terminal time and we will denote by (57, T ^ {Tt)t<^{Vi,T\^ ^ filtered 
probability space, fulfilling the usual conditions. In the whole paper, without restriction of generality T will 
stand for the ct- field Tt- 

2.1 Generating functions 

Let X = {Xt)t(r[o,T] be a real valued stochastic process. 

Definition 2.1. The cumulant generating function of (the law of) Xt is the mapping z i—> Log(E[e^'^']) 
where Log(t(7) = log(|w|) + iArg(w) where Arg(w) is the Argument of w, chosen in ] — tt, tt]; Log is the 
principal value logarithm. In particular we have 

Kxr-D^C with e'^^'f^' = £[6^-^'] , 

where D ;= {z G C | K[e^'^^^'^'^*] < oo, Wt g [0,r]}. In the sequel, when there will be no ambiguity on 
the underlying process X, we will use the shortened notations K.t for W'e observe that D includes the 

imaginary axis. 

Remark 2.2. 1. For all z Cz D, Kt{z) = Kt{z) , where z denotes the conjugate complex of z €z C 
2. For all z e DDR , Kt{z) e M . 
In the whole paper M* will stand for M — {0}. 

2.2 Semimartingales 

An (J't)-semimartingale X = {Xt)teio,T] is a process of the form X ^ M + A, where M is an (J't)-local 
martingale and A is a bounded variation adapted process vanishing at zero. ||^||t will denote the total 
variation of A on [0,T]. If A is (J't)-predictable then X is called an (J't)-special semimartingale. The 
decomposition of an (7't)-special semimartingale is unique, see Definition 4.22 of |25]. Given two {J-'t)- 
locally square integrable martingales M and N, {M, N) will denote the angle bracket of M and N, i.e. the 
unique bounded variation predictable process vanishing at zero such that MN — {M,N) is an (/"t)-local 
martingale. If X and Y are (J't)-semimartingales, [X, F] denotes the square bracket of X and Y, i.e. the 
quadratic covariation of X and Y. In the sequel, if there is no confusion about the underlying filtration (J-'t), 
we will simply speak about semimartingales, special semimartingales, local martingales, martingales. 

All along this paper we will consider C- valued martingales (resp. local martingales, semimartingales). 
Given two C- valued local martingales M^, then M^, are still local martingales. Moreover {M^,AP) = 
(Afi, Af2) . If M is a C- valued martingale then (Af, M) is a real valued increasing process. 

All the local martingales admit a cadlag version. By default, when we speak about local martingales 
we always refer to their cadlag version. Given a real cadlag stochastic process X , the quantity AXt will 
represent the jump Xt — Xt-. More details about previous notions are given in chapter I of |25| . 

For any special semimartingale X we define H-'f ||^2 = E [[M, ATJt] + E (| |^| |t) . The set S'^ is the set of 
(J-t)-special semimartingale X for which ||X||^2 is finite. 
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2.3 Follmer-Schweizer Structure Condition 



Let X = (^t)te[o,T] be a real-valued special semimartiiigale with canonical decomposition, X = AI + A. 
For simplicity, we will just suppose in the sequel that M is a square integrable martingale. For the clarity 
of the reader, we formulate in dimension one, the concepts appearing in the literature, see e.g. |38j in the 
multidimensional case. For a given local martingale M, the space L^{M) consists of all predictable R- valued 
processes v = (wt)tG[o.T] such that E \vs\'^d (M) ^ < oo, where (AI) := {M,AI). For a given predictable 
bounded variation process A, the space L'^{A) consists of all predictable R-valued processes v — {vt)t£[o,T] 
such that E {J^ < oo . Finally, we set 

e := L'^{M)r\L'^{A), (2.1) 

which will be the class of admissible strategies. For any v E Q, the stochastic integral process Gt{v) := 
JgVgdXs, for all t S [0,r] , is therefore well-defined and is a semimartingale in S^. We can view this 
stochastic integral process as the gain process associated with strategy v on the underlying process X. 

The minimization problem we aim to study is the following. Given H € £^(ri,J^, P), a pair {Vq,lp), 
where Vb G M and (p E Q is called optimal if (c, v) — (Vo, (fi) minimizes the expected squared hedging error 

E[{H-c-Gt{v))^] , (2.2) 

over all pairs (c, w) £ R x 0. Vq will represent the initial capital of the hedging portfolio for the contingent 
claim H at time zero. The definition below introduces an important technical condition, see j38) . 

Definition 2.3. Let X = (Xt)tg[o,T] be a real-valued special semimartingale. X is said to satisfy the 
structure condition (SC) if there is a predictable R-valued process a — {at)te[o.T] such that the following 
properties are verified. 

1. At— J^asd{A'I)^ , for all t E [0,T]; in particular dA is absolutely continuous with respect to d{AI) , 
in symbols we denote dA <^ d (AI) . 
T 
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2. / aid{AI)^ < oo , P-a.s. 
Jo 

From now on, we will denote hy K — {Kt)te[o^T] the cadlag process Kt — a^d (M)^ , for all t S [0, T] . 
This process will be called the mean-variance trade-off (MVT) process. Lemma 2 of |38| states the 
following. 

Proposition 2.4. If X satisfies (SC) such that Kt is a bounded r.v., then Q = L^(Af). 

The structure condition (SC) appears naturally in applications to financial mathematics. In fact, it is 
mildly related to the no arbitrage condition at least when X is a continuous process. Indeed, in the case 
where X is a continuous martingale under an equivalent probability measure, then (SC) is fulfilled. 



2.4 Follmer-Schweizer Decomposition and variance optimal hedging 

Throughout this section, as in Section [^751 X is supposed to be an (J^t)-special semimartingale fulfilling the 
(SC) condition. Two (J^t)-martingales AI,N are said to be strongly orthogonal if AdN is a uniformly 
integrable martingale, see Chapter IV. 3 p. 179 of [31]. If M,N are two square integrable martingales, then 
M and iV are strongly orthogonal if and only if {AI, N) — 0. This can be proved using Lemma IV. 3. 2 of |31| . 
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Definition 2.5. A random variable H G £^(r2, J^,P) admits a Fdllmer-Schweizer (FS) decomposition, 
if 

H = Hq+ f ^^dX, + , P~ a.s. , (2.3) 
Jo 

where Hq G M. is a constant, E Q and = (L^)(£[o,t] square integrahle martingale, with W^L^] = 
and strongly orthogonal to M . 

We summarize now some fundamental results stated in Theorems 3.4 and 4.6, of |30| on the existence 
and uniqueness of the FS decomposition and of solutions for the optimization problem 



Theorem 2.6. We suppose that X satisfies (SC) and that the MVT process K is uniformly hounded in t 
andoj. Let H e C'^{n,T,P). 

1. H admits a FS decomposition. It is unique in the sense that Hq G M, £ L'^{M) and is uniquely 
determined by H . 

2. For every H e C^{n,F,V) and every c G £^(J"o), th ere exists a unique strategy ip^'^'^'^ £ O .such that 

E[{H - c - GT{<p^^'"^)f] = minE[(i7 - c - Gt(iO)'] • (2-4) 

3. For every H G £^(fi, J-", P) there exists a unique couple {c^^\lp^^'>) G £'^(J^o) x 9 such that 

E[(i7-c(«)_GT(^(^'))^]-, nrm mH - c - Gt{v)Y] . 

Next theorem gives the explicit form of the optimal strategy under some restrictions on K. 

Theorem 2.7. Suppose that X satisfies (SC) and that the MVT process K of X is deterministic and let a 
he the process appearing in DehnitionlEE of (SC). Let H G C^{n,F,V) with FS decomposition ifO)) . 

1. For any c G M, the solution of the minimization problem {2.4^ verifies Lp^'^'^^ G Q, such that 

vf'"^ - ef + YrkK^'- -^-Gt^iv^'-"^)) ' V all t G [0,T] (2.5) 
where the process {Ht)te[o.T] *s defined by Ht := Hq + j^^^^dXg + . 

2. The solution of the minimization problem 12. 2\) is given by the pair (i/o, <y3^^°'^'') ■ 

3. If (M) is continuous, 
imE[{H~c~GT{v)f] = cxp{~KT){iHo-cf +E[{L^f]) 



mml 



+E 



r cM-{KT-K.,)}d{L")^ 
Jo 



Proof. Item 1. is stated in Theorem 3 of [S^- Item 2. is a consequence of Corollary 10 of [38]. Item 3. is a 
consequence of Corollary 9 of [3H] taking into account that K inherits the continuity property of (M). We 
remark that K = K, where iiT is a process appearing in the statement of the mentioned corollary. □ 

In the sequel, we will find an explicit expression of the FS decomposition for a large class of square 
integrable random variables, when the underlying process is an exponential of additive process. 
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2.5 Additive processes 

This subsection deals with processes with independent increments which are continuous in probabihty. From 
now on (J-'t) will always be the canonical filtration associated with X . 

Definition 2.8. A cadlag process X ~ (-''^t)tG[o.T] is a (real) additive process iff Xq — 0, X is continuous 
in probability, i.e. X has no fixed time of discontinuities and it has independent increments in the following 
sense: Xt — Xs is independent of Ts for < s < t <T. 

X is called Levy process if it is additive and the distribution of Xt — Xg only depends on t — s for 
0<s<t<T. 

An important notion, in the theory of semimartingales, is the notion of characteristics, introduced in 
definition II. 2. 6 of [53]. A triplet of characteristics (6, c, depends on a fixed truncation function h : 
M — > K. with compact support such that h{x) = a; in a neighborhood of 0; is some random cr-finite Borel 
measure on [0, T] x M. If X is a semimartingale additive process the triplet (6, c, v) admits a deterministic 
version, see Theorem II. 4. 15 of Moreover (bt), (q) and t i— > /[otjxsd^l^ ^ l)i^{ds,dx) have bounded 
variation for any Borel real subset B. Generally in this paper B{E) denotes the Borel cr-field associated with 
a topological space E. 

Proposition 2.9. Suppose X is a semimartingale additive process with characteristics (6, c, i/), where v is 
a non-negative Borel measure on [0,T] x R. Then t^ at given by 

at = ||6||t + Q+ I {\x\^ M)v{\Q,tldx) (2.6) 

fulfills 

dbt dat , dct <^ dat and v{dt^dx) = Ft{dx)dat , (2-7) 
where Ft{dx) is a non-negative kernel from ([0, T], S([0, T])) into (R, i3(M)) verifying 

I {\x\^ Al) Ft{dx) <1, V<e[0,T]. (2.8) 

JR 

Proof. The existence of (at) as a process fulfilling (|2.6[) and F fulfilling (|2.8p is provided by the statement 
and the proof of Proposition II. 2.9 of [IS]. (|2.6p guarantees that (at) is deterministic. □ 

We come back to the cumulant generating function k and its domain D. 

Remark 2.10. In the case where the underlying process X is an additive process, then 

D:={zeC\ E[e-^"(^)-^'] < oo, Vi G [0,T]} = {z e C \ E[e'^''^'^^^] < oo} . 

In fact, for given t £ [0,r],7 e M we have E(e''-^'^) = E(e''-^')E(e'^('^'^"'^'')) < oo. Since each factor is 
positive, if the left-hand side is finite, then E(e''' '''"') is also finite. 

3 FoUmer-Schweizer decomposition for exponential of additive pro- 
cesses 

The aim of this section is to derive a quasi-explicit formula of the FS decomposition for exponential of 
additive processes with possibly non stationary increments. 



7 



We assume that the process S is the discounted price of the non-dividend paying stock which is supposed 
to be of the form, St = SQCxjp{Xt), for all t € [0,T], where sq is a strictly positive constant and X is a 
semimartingale additive process, in the sense of Definition l2.81 but not necessarily with stationary increments. 
In the whole paper, if z is a complex number, stands for exp(ln(so) + zXt). In particular if y is a real 
number, stands for soexp{yXt). 

3.1 On some properties of cumulant generating functions 

We need now a result which extends the classical Levy-Khinchine decomposition, see e.g. 2.1 in Chapter II 
and Theorem 4.15 of Chapter II, which is only defined in the imaginary axis to the whole domain of the 
cumulant generating function. Similarly to Theorem 25.17 of )35| . applicable for the Levy case, for additive 
processes we have the following. 

Proposition 3.1. Let X be a semimartingale additive process and set Dq = Cz R \ J[q x]x{\x\>i} ^'^'^'^^^^^'^^') < oo 
Then, 

1. Dq is convex and contains the origin. 

2. Do^ DnR. 

3. If z Cz C such that Re{z) € Dq, i.e. z £ D, then 

z^ f 

Kt{z) = zbt + —ct+ / (6^^= - 1 - zh{x))i^{ds, dx) . (3.1) 

Proof. 1. is a consequence of Holder inequality similarly as i) in Theorem 25.17 of |35j . 

2. The characteristic function of the law of Xt is given through the characteristics of X , i.e. 

*t(u) = iubt ct + / (e™^ - 1 - iuh{x))F\dx) , for aU u G R, 

2 JR 

where we recall that for any t > 0, q > and B H> F*{B) := v{[Q, t] x i?) is a positive measure which 
integrates 1 A |xp. Let t G [0,T]. According to Theorem II. 8.1 (iii) of [3S] , there is an infinitely divisible 
distribution with characteristics ct,F*{dx)). By uniqueness of the characteristic function, that law 
is precisely the law of Xt. By Corollary II. 11. 6, in [35]; there is a Levy process (^^,0 < s < 1) such 
that L{ and Xt are identically distributed. We define 

C^ = {ceM| / e'''' Ft{dx) <oo} and C* {z e C | E [cxp(i?e(zL*i)] < oo} . 

J{\x\>l} 

Remark 12.101 says that = D, moreover clearly = Dq. Theorem V.25.17 of [35] imphes Dq = 
D DM., i.e. point 2. is established. 

3. Let t G [0,r] be fixed; let z E D C C*, in particular Re{z) G Cg. We apply point (iii) of Theorem 
V.25.17 of [35] to the Levy process L*. 

□ 

Proposition 3.2. Let X be a semimartingale additive process. For all z E D, t t-^ Ht{z) has bounded 
variation and K.dt{z) ^ rfa*, where t t-^ at was defined in Provosition \2.{A 
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Proof. Using p.ip . we only have to prove that t H> t] xr(^^^ — 1 — zh{x))h'{ds, dx) is absolutely continuous 
w.r.t. {dat). We can conclude 

/■* dbs /■* dcs /■* [ 

Kt{z) = z I -f^das + Y / T^das + dos (e^^ - 1 - zh{x)) Fs{dx) , 



if we show that 



das / le^"^ - 1 - z/i(a;)|F,(da;) < oo . (3.2) 

Without restriction of generality we can suppose h{x) ~ xl^^^^^i. (j3.2p can be bounded by the sum I1+I2 + I3 
where 

Ii= [ das [ |e""|i^,(da;) , h = [ das [ Fs{dx) , and I3 = [ das ( \e'^'=-l-zx\Fs{dx) . 

Jo ^Ix^l Jo J\x\>l Jo "'|a;|<l 

Using Proposition 12.91 we have 

Ii= [ das [ \e''=\Fs{dx)= [ das ( e^'^^^^^F^dx) = /" e^'=(^)^i^(ds, dx); 

Jo ^|3;|>l Jo J\x\>l J[0,T]x{|2:|>l} 

this quantity is finite because Re{z) G Dq taking into account Proposition 13. II Concerning I2 we have 
h= I das ( Fs{dx)= [ das j {I h\x^\)Fs{dx) < ar, 

Jo "'|x|>l Jo J\x\>l 

because of (|2.8p . As far as is concerned, we have 

h < e«^(^)^ / das{x^ A l)Fs{dx) ^ e«^(^)^ar 

^ J[0,T]x{|x|<l} ^ 

again because of (|2.8p . This concludes the proof of the proposition. 

□ 



The converse of the first part of previous Proposition 13.21 also holds. To show this, we formulate first a 
simple remark. 

Remark 3.3. 1. For every z Cz D, {exjp{zXt — Kt{z))) is a martingale. In fact, for all < s < t < T , we 
have E[exp(z(Xt — ^s))] = cxp(Kt(z) — Ks{z)). 

2. t 1-^ Kt{0) = 1 and it has always bounded variation. 

Proposition 3.4. Let X be an additive process and z £ ZJnM*. {Xt)ti=[o^T] ci semimartingale if and only 
if t 1-^ Kt{z) has bounded variation. 

Proof. Using Proposition 13.21 it remains to prove the converse implication. If i 1— > Kt{z) has bounded 
variation then t 1— > e'^'^^' has the same property. Remark 13.31 savs that e^-^* = A/te"'*^^^ where (Mt) is a 
martingale. Finally, {e'^'^*) is a semimartingale and taking the logarithm (zXt) has the same property. □ 

Remark 3.5. Let z G D. If (Xt) is a semimartingale additive process, then (e^-^*) is necessarily a special 
semimartingale since it is the product of a martingale and a bounded variation continuous deterministic 
function and by use of integration by parts. 

Proposition 3.6. The function {t,z) 1— >■ K,t{z) is continuous. In particular, {t,z) 1— >■ Kt{z), t € [0,2"], z 
belonging to a compact real subset, is bounded. 
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Proof. • Proposition 13. II implies that z i— >■ Kt(z) is continuous uniformly w.r.t. t G [0,T]. 

• We first prove that z E lnt{D), t i-> Kt{z) is continuous. Since z G fnt(£'), there is 7 > f such that 
72 G D; so 

E[exp{zjXt)] = exp(Kt(7z)) < exp(sup(Kt(7z))) , 

t<T 

because t 1— > Kt{'jz) is bounded, being of bounded variation. This imphes that (exp(zXi))jg[o,T] is 
uniformly integrable. Since (X^) is continuous in probability, then {exp{zXt)) is continuous in C^. The 
partial result easily follows. 

• To conclude it remains to show that 1 1— )■ Kt{z) is continuous for every z Cz D. Since D = Int(£'), there 
is a sequence (z„) in the interior of D converging to z. Since a uniform limit of continuous functions 
on [0, T] is a continuous function, the result follows. 

□ 

3.2 A reference variance measure 

For notational convenience we introduce the set ^ = {z G C\2z G D}. 
Remark 3.7. We recall that D is convex. Consequently we have. 

1. Ify,ze -J, then y + zeD. If z e ^ ^^e*^ ^ G T "'^'^ 2Re{z) G D. 

2. Since e D, clearly C £>. 

3. Under Assumption[J\ below, 2 G -D and so + 1 C -D. 

We introduce a new function that will be useful in the sequel. 
Definition 3.8. • For any t G [0,T], if z,y £ -j we denote 

pt{z, y) = Kt{z + y)- Kt(z) - Kt{y) . (3.3) 

• To shorten notations pt : ^ ^ C will denote the real valued function such that, 

pt{z) = pt{z, z) = Kt{2Re{z)) - 2Re{Kt[,z)) . (3.4) 

Notice that the latter equality results from Remark \2.2\ 1. 

An important technical lemma follows below. 

Lemma 3.9. Let z G -j, with Re{z) ^ 0, then, t H- pt{z) is strictly increasing if and only if X has no 
deterministic increments. 

Proof. It is enough to show that X has no deterministic increments if and only if for any < s < t < T, the 
following quantity is positive, 

Pt{z) - ps{z) = [Kt{2Re{z)) - K,{2Re{z))] - 2Re{nt{z) - At,(z)) . (3.5) 

By Remark l3.3i we have eyiY)[Kt{z) — Ks{z)] = E[cxp(zA*X)], where A*X := Xt — Xg. Applying this property 
and Remark 12.21 1., to the exponential of the first term on the right-hand side of p.Sp yields 

exp [Kt{2Re{z)) - K,{2Re{z))] = E[exp(2i?e(z)A*X)] = E[cxp((z + z)A*X)] = E[|exp(zA*X)|^] . 
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Similarly, for the exponential of the second term on the right-hand side of p.Sp . one gets 

exp [2Re(Kt{z) - Ksiz))] = cxp (Kt{z) - Ks{z)) + (Kt{z) - Ks(z)) = |E[exp(zA*X)]|^ 
Hence taking the exponential of A*p(z) := pt{z) — Ps{z) yields 



TAt -, E[|exp(zA*X)|^] E[|r*X(z)|"] ^ , ^ , w 

exp A*p z) - 1 = _LL_tli " ^' ^ - 1 = -LL^ — - l , where T\X{z) = exp(zA*X) 

'E[exp(zA*X)]|2 m\X{z)f 

(3.6) 

"t^ar [i?e(r*X(z))] + Var [lm{V\X(z))\ 

\nnx{z)]\^ ■ 

• If X has a deterministic increment A*X = Xt — Xs^ then r*X(2:) is again deterministic and p.6p 
vanishes and hence < H> pt(z) is not strictly increasing. 

• li X has never deterministic increments, then the nominator is never zero, otherwise Re{T\X{z)^ = 
exp(i?e(z)A*X), and therefore A*X would be deterministic. 

□ 

^ ^ , ^ l/ar(exp(A*X)) 

Remark 3.10. 1/2 G 13, setting z ~ \ in hS. b]) implies that ptyVj > Ps(l) is equivalent to 5- > 

(E[exp(AiX)]) 

0. Taking the process S at discrete instants to ~ < ■ ■ ■ < tk < ■ ■ ■ < tn ~ T , one can define the discrete 
time process {S^)k=o,... such that S'^ ~ Stj. and derive the counterpart of Lemma \S.9\ in the discrete time 
setting. Indeed, the following assertions are equivalent: 

• (Ptfc (l))fc=o,---n o,n increasing sequence; 

• Aj^^'^X is never deterministic for any k = 0, ■ ■ ■ ,n — 1. 

Moreover, accordingly to Proposition 3.10 in ]2'^ . we observe that, under one of the above equivalent condi- 
tions, the (discrete time) mean-variance trade-off process associated with (5'fc)fc=o.- - ,n defined by 

^ g {nSk+i'Su\Tu]f ^ ^ g (E[exp(A:r^^) - if 



{Var[Sk+i - Sk\Tk]) {Var [exp(A*^+^X)J f 

is always bounded. According to Proposition 2.6 of l^Oj , that condition guarantees that every square integrable 
random variable admits a discrete Fdllmer-Schweizer decomposition. The process K'^ is the discrete analogous 
of the MVT process K; one can compare the mentioned result to item 1. of Theorem \2.b\ 

From now on, we will always suppose the following assumption. 

Assumption 1. 1. {Xt) has no deterministic increments. 

2. 2eD. 

We continue with a simple observation. 
Lemma 3.11. Let I be a compact real interval included in D. Then sup^.^/ sup(<2i E[S'f ] < 00. 
Proof. Let i € [0, T] and x Cz I, since k is continuous, we have 

E[S'f] = sgexp{Kt(a;)} < max(l, 4"^^) exp(supt<y^^gj \nt{x)\) . □ 
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Remark 3.12. From now on, in this section, dpt = pdt will denote the measure 



dpt^ Pdt{l)^d{Kt{2)-2nt{l)) . 



(3.7) 



According to Assumptions^ and Lemma \3.9\ it is a positive measure which is strictly positive on each interval. 
This measure will play a fundamental role. 

We state below a result that will help us to show that K.dt{z) is absolutely continuous w.r.t. pdti^)- 

Lemma 3.13. We consider two positive finite non-atomic Borel measures on E (Z M", p and v. We suppose 
the following: 

1. p <^ V ; 

2. p{I) 7^ for every open ball I of E. 

Then h := — ^ ^ Q v a.e. In particular p and v are equivalent, 
dv 

Proof. We consider the Borel set B = {x ^ E\h{x) = 0}. We want to prove that v{B) = 0. So we suppose 
that there exists a constant c > such that i^^B) ~ c > and take another constant e such that < e < c. 
Since v is a Radon measure, there are compact subsets K,, and of E such that C K± C B and 
i'{B — Kf) < e, v{B — K^) < |. Setting e = |, we have v{Kf) > | and v{K^) > ^. By Urysohn lemma, 
there is a continuous function (/s : _B — > R such that, < < 1 with = 1 on and = on the closure 
of K%. Now if{x)v{dx) > v{K^) > | > 0. By continuity of f there is an open set O C E' with 1^9(2;) > 
for a; G O. Clearly O C K^. C B; since O is relatively compact, it is a countable union of balls, and so 
B contains a ball /. The fact that h = Q on I implies p{I) = and this contradicts Hypothesis 2. of the 
statement. Hence the result follows. 



Remark 3.14. 1. If E = [0,r], then point 2. of Lemma \3.13\ becomes p{I) 7^ for every open interval 



2. The result holds for every normal metric locally connected space E, provided v are Radon measures. 
Proposition 3.15. Under Assumptionl^ 



Proof. We apply Lemma 13.131 with dp = dpt and dv = dat. Indeed, Proposition 13.21 implies Condition 1. of 
Lemma 13.131 and Lemma 13.91 implies Condition 2. of Lemma 13. 131 Therefore, dat is equivalent to dpt . □ 



□ 



/ c [0, r 



d(Kt{z)) <^ dpt , for all z Cz D . 



(3.8) 



Remark 3.16. Notice that this result also holds with dpt{y) instead of dpt = dpt{l), for any y G 
that Re{y) ^ 0. 



■J such 



3.3 On some semimartingale decompositions and covariations 

Proposition 3.17. We suppose the validity of item 2. of Assumption{l\ Lety,z G -j. Then is a special 
semimartingale whose canonical decomposition = M{z)t -\- A[z)t satisfies 




(3.9) 



where dpu(z) is defined by equation ^3^. In particular we have the following: 
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1. {M{z),M), = f^Slt^pau{z,l) 

Remark 3.18. • Clearly 1 e Z? since and 2 belong to Dq and Dq is convex by Provosition \S.l\ 

• If z = 1, we have = S, so that by uniqueness of the special semimartingale decomposition, it follows 
that A/(l) = M. 

Proof. The case y = 1, follows very similarly to the proof of Lemma 3.2 of [231 • The major tools are 
integration by parts and Remark 13.31 which says that N{z)t :— e~'''(^^5f is a martingale. The general case 
can be easily adapted. □ 



Remark 3.19. Lemma \3.11\ imvlies that E [| (M(y), Af(z)) |] < oo and so M{z) is a square integrable mar- 
D 

2 



tingale for any z G ^ 



3.4 On the Structure Condition 

Proposition 13.171 with y = z = 1 yields S M + A where At = S u—t^du{^) and AI is a martingale such 
that {M,M)^ ~ Jo S'u^i^duC^) — 2K(i„(l)) = Jq S^-Pdu- At this point, the aim is to exhibit a predictable 
R- valued process a such that 

1. At = J*asd{M)^,te[0,T]. 

2. Kt = ald{M)^ is bounded. 

In that case, according to item 1. of Theorem 12. 6[ there will exist a unique FS decomposition for any H G 
£^(r2, J^, P) and so the minimization problem (|2.2p will have a unique solution, characterized by Theorem 12. 71 
2. 

Proposition 3.20. Under Assumptions^ At = agd {M) ^ , where a is given by 

" - ^" with A„:=%^, for allu(,[Q,T]. (3.10) 



bu- dp. 
Moreover the MVT process is given by 



Kt= I (3.11) 



\ 'iPu 

Corollary 3.21. Under Assumption\^ the structure condition (SC) is verified if and only if 



Jo \ dpu J 
In particular, (Kt) is deterministic therefore bounded. 

Remark 3.22. Item 1. of Assumptions^ is natural. Indeed if it were not realized, i.e. if X admits a 
deterministic increment on some interval [s,t], then S would not fulfill the (SC) condition, unless u i— > k„(1) 
is constant on [s,t\. In this case, the market model would admit arbitrage opportunities. 

Proof (of Proposition [3. 20p . By Proposition 13. 15) dKt{l) is absolutely continuous w.r.t. dpt. Setting as 
in (|3.10p . relation p. lip follows from Proposition 13. 171 expressing Kf = j^a\d{M)^. □ 
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Lemma 3.23. The space Q, defined in (|2.ip . is constituted by all predictable processes v such thatW. ^ Vf Sf_dpt^^ 



< 



00. 



Proof. According to Proposition [2111 the fact that K is bounded and S satisfies (SC), then v E Q holds if and 



r 


dKu{z) 


Jo 


dpu 



only if V is predictable and E[Jq v^d{M,M)^] < oo. Since {M,M)^ ~ S'^_dps, the assertion follows. □ 

3.5 Explicit Follmer-Schweizer decomposition 

We denote by V the set oi z E D such that 

dpu < oo. (3-12) 



From now on, we formulate another assumption. 
Assumption 2. lei?. 

Remark 3.24. 1. Because of Provosition W.15[ — exists for every z E D. 

apt 

2. Under Assumption\^ Corollarv \3.21\ says that Assumption\^is equivalent to (SC). 

The proposition below will constitute an important step for determining the FS decomposition of the 
contingent claim H ~ fiSr ) for a significant class of functions /, see Section [5751 

Proposition 3.25. Let z G 2? n with z + 1 £1), (in particular 2Re{z) £ D), then 

2. Moreover, suppose that Assumptions\^and\^hold and define 

d{pt{z,l)) 



dpt 



te[o,T] 



Then \^iz,t)\'^p^f < oo and 



/"* dKs(l) 

r]{z,t) := Kt{z) - I 7(z,s)Kds(l) = Kt(2) - / ^{z,s) 



Pds 



(3.13) 



(3.14) 



/ Jo dps 

is well-defined, besides ri{z,-) is absolutely continuous w.r.t. pds and therefore bounded. 

3. Again under Assumptions\^and\^ H(z) ~ admits an FS decomposition H{z) = H{z)o+ ^{z)tdSt + 
L{z)t where 



(3.15) 
(3.16) 

(3.17) 



az)t 7(^,06^""^^^'^^^*-"', 

L{z)t H{z)t - i/(z)o - / az)udSu 

Jo 

Proof. 1 . is a consequence of Lemma 13.111 

2. 7(2;, •) is square integrable because Assumption [2l and z, z + 1 G P. Moreover rj is well-defined since 

^^^(1) 



h{z,s)\ 



dps 



Pds] < \liz,s)\ Pds 



f 


dKs{l) 
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dps 



Pds- 



(3.18) 
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3. In order to prove that (|3.15p . (|3.16p and (|3.17p is the FS decomposition of H{z), we need to show that 

(a) H{z)q is J^o-iT^iea-surable, 

(b) (i(z),A/) =0, 

(c) £,{z) e e, where 6 was defined in (|2.ip . 

(d) L{z) is a square integrable martingale. 

We proceed similarly to the proof of Lemma 3.3 of fMj. Point (a) is obvious. Partial integration and 
point [T] of Proposition 13.171 yield 



H{z)t^H{z)o+ eLv{^^'i^)dM{z)u- ^i^^^^) Sli^{z,du) + eJ. vi.^^'^-) Sl_Kdu{z) . (3.19) 

Jo Jo Jo 

On the other hand 

/ £,{z)udSu = I az)udMu+ f j{z,u)ef^'^^'^'''^S:,_Kdu{l) ■ (3.20) 
Jo Jo Jo 

Hence, using expressions p.lPp and p.20p . by definition of rj in (j3.14p . which says 77(2, du) — Hduiz) — 
7(2, u)«;d„(l), we obtain 

L{z)t = H{z)t-H{z)o- [ i{z)udSu= I e^" "(^''*^)dAf(z)„ - / (3.21) 



which implies that L{z) is a local martingale. 

From point [TJ of Proposition I3.17[ using (|3.16p , it follows that 

(L(z),M), = f e^^'^^'^'''^Slt\pdu{zA)--i{z,u)pdu] 



Then by definition of 7 in (j3.13l) . pdt{z, 1) — "f{z,t)pdt , yields {L(z), M)^ = 0. Consequently, point (b) 
follows. 

It remains to prove point (d) i.e. that L(z) is a square-integrable martingale for all z G -D and that 
i?e(^(z)) and Im{£^{z)) are in 8. (|3.2ip says that 

L{z)t^ re/"''(^-*")dAf,(z)- f\{z)sdM,. 
Jo Jo 

By Remark |2. 21 we observe first that z, z + 1 G P. Moreover by definition of 7 and rj. it follows 



7(z,t) = 7(z,t) and ri{z,t) ^ ■qiz.t). (3.22) 
By Proposition [STfi WT2\ and it follows 

(l{z),L{z))^ - (L(z),L(z-)),==^L(z),^'e/""(^"-*")dM,(z)^ 



t 

t 



Consequently 

(L(z),IR)^=^^/=^^^■^(''(-^"))5,^^^W[p..(z)-|7 (3.24) 
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Taking the expectation in p.24p . using point 2., (|3.13p . (|3.14p and Lemma [3.111 we obtain 

E 



Liz), Liz)) <oo. (3.25) 
/ tJ 

Therefore, L is a square-integrable martingale. 

It remains to prove point (c) i.e. that ^(z) € 6. In view of applying Lemma 13.231 we evaluate 



JO 



Similarly as for (|3.24p , we can show that the expectation of the right-hand side of p.26p is finite. This 
concludes the proof of Proposition 13.251 

□ 

3.6 FS decomposition of special contingent claims 

We consider now payoff functions of the type 

H = fiSr) with /(,s) - / s'Uidz) , (3.27) 

JC 

where 11 is a (finite) complex measure in the sense of Rudin |33) . Section 6.1. An integral representation of 
some basic European calls is provided in the sequel. We need now the new following assumption. 

Assumption 3. Let Iq = suppll n M. We denote I — 2Iq U {1}. 

L Iq is compact. 

2. Vz e suppII, z,z + leV. 



3. h^i 



4- sup^gj 



d{Kt{x)) 
dpt 



< oo. 



Remark 3.26. L Two kinds of assumptions appear. Assumptions]^ and\^ only concern the process and 
Assumption]^ involves both the process and the payoff. 

2. Assumption\3l looks obscure. Examples for its validity will be provided in Section]^ For instance con- 
sider the specific case where X is a Wiener integral driven by a Levy process A, i.e. Xt = /g l{s)dAs, t e 
[0, T] and the payoffs are either a call or a put. We observe in Example \5.6\ below that Assumptions]^]^ 
and\^are a consequence of the simple Assumption^ 

Remark 3.27. L Point 3. of Assumption]^ implies sup^g^^jjj ||K£it(i?e(z))||y < oo . 

2. Under Assumption\^ H = /(S't) is square integrable. In particular it admits an FS decomposition. 

3. Because of 113. 8\} in Provosition \3.15[ the Radon-Nykodim derivative at Assumption]3i4, always exists. 

We need now to obtain upper bounds on z for the quantity (|3.25p . We will first need the following lemma 
which constitutes a (not straightforward) generalization of Lemma 3.4 of |24| which was stated when X is a 
Levy process. The fact that X does not have stationary increments, constitutes a significant obstacle. 

Lemma 3.28. Under Assumptions]^]^]^ there are positive constants Ci,C2,C3 such that dps a.e. 
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dReifilz, s)) 

1. sup — y < ci . 

2. For any z G /q + , |7(z, s)p < '^^^^^^ 



< C2 - C3 



dRe{r/{z, s)) 



dps 



3. 



sup / 2Re{T]{z, dt)) cxp ( [ 2Re{7]{z,ds))] < 

zelo+iRJo \Jt J 



Remark 3.29. 1. According to Provosition \3.25[ t Re{ri{z,t)) is absolutely continuous w.r.t. dpt 
2. We recall that suppll is included in Iq + iR. 

Proof (of Lemma 13. 28p . According to Point 3. of Assumption [3] we denote 

d{Kt{x)) 



Oil := sup 



dpt 



(3.28) 



For z G /q + iM, t e [0, T], we have 77(2, t) = nt{z) — 7(2, s)o?Ks(1) and 77(2, t) = Kt{z) — 7(2, s)dKsi^) ■ 
Then, we get Re{r]{z,t)) = Re{Kt{z)) — Re{'-f{z, s))dKs{l) ■ We obtain 



Re{T]{z,ds)) < Re{KT{z) - Kt{z)) + 
'■^ Re{dKs{z)) 



dps 



-dps 



j{z,s)dKs{l) 



7(z,s)dKs(l) 



(3.29) 



Since (^L{z),L{z)j is increasing, taking into account p.24p . the measure (dps{z) — \-^{z, s)\'^dps) is non- 
negative. It follows that 



dpsjz) 
dps 



\l{z, s)p > , dps a.e. 



(3.30) 



By (|3.30p . in particular the density ^ jg non- negative dps a.e. Consequently, 

dps 

dRe{K,(z)) dKsi2Re{z)) 

2 ; < , dps a.e. 

dps dps 

In order to prove 1. it is enough to verify that, for some co > 0, 

dReirjiz^s)) 1 dRe(Ks(z)) 

< Co + ; dps a.e. 



dps 



2 dps 



(3.31) 



(3.32) 



dRe(v(z s)) 1 

In fact, (|3?3T]) . Assumption[3]point 3. and (|3?28l) . imply that ' < ca + -cn =: ci. To prove (|332l) 

dps 2 

it is enough to show that 



Re{Tj{z,T) - 7j{z,t)) < coipT ~ Pt) + ^Re{KT{z) - Kt{z)), \ft e [0,T]. 



Again Assumption [3] point 3. implies that 



'-f{z,s)dKs{l) 



< ci: 



• j^ \liz,s)\dps , 



(3.33) 



(3.34) 
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where C12 



Using p.30[) and Assumption [3] it follows 



, , ,,2 / dp,{z) dK.{2Re{z)) 2dRe{Ks{z)) ^ 2dRe{K,{z)) 

h{z,S)\ <— = ; ; < Cll 



dps 



dps 



dps 



dps 



(3.35) 



This implies that c^2 l7(-2,s)|^ < 
0^2 Cll- Consequently, 



(c?3 + \ ( 



where C13 > is chosen such that cfg > 4C]^2 



7(2,5)^^5(1) < / ( Ci3 

Coming back to p.29p . we obtain 

i?e(?7(z,T)-77(z,t)) < 

< 

(|3.3ip and Assumption [3] allow to establish 



dRe{Ks{z)) 



dps 



f Re{dKs{z)) 1 



Re{dKs{z)) 



dps 



1 Re{dhis{z)) ( Re{dKs{z)) 



dps 



dps 



+ ci3 



R r r ^ IdReMf)) 
Re{ri{z,T)~r){z,t)) < dps Icq + ^ ^ 

where cq = ^ + C13. This concludes the proof of point 1. 
In order to prove point 2. we first observe that (|3.32p implies 



(3.36) 



dRe{Ks{z)) 



dps 



< 2 Co 



dRe{ri{z, s)) 
dps 



dps a.e. 



(3.37) 



p.35p implies \j{z, s)|^ < C21 — 4 



dRe{r]{z, s)) 
dps 



where C21 = cn + Icq. Point 2. is now established with 



C2 = C21 and C3 = 4. 

We continue with the proof of point 3. We decompose Re{r]{z,t)) = A'^{z,t) — A^{z,t), where A^{z, .) and 
A~(z, .) are the increasing non negative functions given by 



A+{z,t) 



( dRe{r]{z, s)) 



\ dps 



(dRe{r^{z,s)) 



V dps 



dps 



dps and A (z, t) 

Moreover point 1. implies j4+(z,i) < cipt- At this point, for z e /q + iM 

i?e(77(z, dt))ef" 2i?e(r,(.,ds)) ^ - ^+ (z, d^) R^M^^^s)) 



< 



< e 



2ci/: 



g2cipr 



which concludes the proof of point 3 of Lemma 13.281 



|l_e-2A-(..T)| < 



□ 



Theorem 3.30. Let IT he a finite complex- valued Borel measure on C. Suppose Assumptions\^\^\^ Any 
complex-valued contingent claim H = /{St), where f is of the form i3.27\) . and H G £^(51, J^, P), admits a 
unique FS decomposition H = Hq + dSt + Lf^ with the following properties. 
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1. H e C^in^T, 



Ht^ I H{z)tll{dz), = J e(^)tn(dz), if = y L{z)tllidz), 



where for z S supp(n), H{z),S,{z) and L{z) are the same as those introduced in Provosition \8. 25\ and 
we convene that they vanish if z ^ supp(n). 

2. Previous decomposition is real-valued if f is real-valued. 
Remark 3.31. Taking 11 = Sz„{dz), zq e C, Assumption\^ is equivalent to the assumptions of Proposition 

Proof (ofTheorem[3l30l). a) fiSr) G C^{n,T,F) since by Jensen's, E \J^Il(dz)S^\'^ < \n\{dz)E\Sl^''^"^\\Il\{C) < 
sup^^jg £'(S'|?)|n|(C)^, where |n| denotes the total variation of the finite measure 11. Previous quantity is 
bounded because of Lemma 13.111 

b) We go on with the FS decomposition. We would like to prove first that H and are well defined 
square-integrable processes and E{J^ \^^\'^d{M)s) < oo. 
By Jensen's inequality, we have 



E 



< 



L{z)tn{dz 
Similar calculations allow to show that 

E[(e^)?] < |n|(C) / \U\dz)Emz)t 
Jc 

We will show now that 

• (Al): supj<j,^,g3,ppnE[|i?t(^)n < 00 
. (A2): /jn|(dz)E[|L(z)Tp]<oo; 



E ^ mdz)\L{z)tf'^ |n(c)| = y |n|(dz)E[|L(z),|2]|n(c)|. 



and Epf )^] < |n(C)| / \Il\idz)E[\Liz)t\ 



. (A3): dp,SU^Uz)midz)) <oo. 

(Al): Since H{z)t = e^T vi^^ds) gz ^ ^^^^ \H{z)t\^ = H{z)t'H{z)'t = e^t"' 2fle(„(^,ds))52i?e(z) 

E[|i7(z)t|2] = 2Re{r,{.As))^^g2Re(z)^ ^ ^^^J^ 2Re(r,(z,ds)) ^ 

where C4 is well defined by p.38p . below, since by Lemma [3. Ill 

C4 := sup E [S"^] < 00 . 

x£l,s<T 



(3.38) 



Lemma [3.281 implies (Al). Therefore [Ht) is a well-defined square-integrable process. (A2): E[|Lt(2)p] < 
E[|Lt(^)P] = E[(L(z), L{z)) ] , where the first inequality is due to the fact that |Lt(z)p is a submartingale. 



E 



Hz), Liz) 



= E 



^f-2R.MzMu))s2Reiz) ^^^^^^^ _ |^(,^,)|2rf^^ 



By Fubini's theorem, Lemma [3.111 and p.24p . we have 



L{z),L{z) 



Sj 2Re(i-,(z,du))^ ^g2Re(z) 



< C4 



,/j2i?e(r,(z,du)) 



dps 



dpsjz) 
dps 

dps- 



hiz,s)\' 



dps 
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According to Lemma 13.281 point 2, previous expression is bounded by 04/(2), where 



' ''^ dRe{T]{z,t)) 



dpt exp / 2Re{7]{z , ds)) 



C2 - C3- 



dpt 



(3.39) 



where Ii{z) = dpt exp (^J^ 2Re{ri{z , ds))j and hiz) — — cxp (^J^ 2Re{ri{z, ds))j Re{ri{z,ds)). Using 
again Lemma [3.281 we obtain 



sup |/i(z)| < pt exp {2ciPt) and sup 1/2(2^)! < 00 , 



and so 



sup 



Liz), Liz) 



< 00 



(3.40) 



(3.41) 



This concludes (A2). 

We verify now the validity of (A3). This requires to control 



E 



l^.tS? (^ljU\idz)\az)t\'^ <E j\dtSUj^mdz) ^iz,t)expij^ 



7(z,t)exp I / i?e(?7(z,ds)) I SI 



Using Jensen's inequality, this is smaller or equal than 



n(c)| / \n\idz) / 



s. 



2Re(z) 



|7(2,t)pcxp 2 / i?e(r/(z,ds)) 



Lemma (3.281 gives the upper bound C4|n|(C) |n|((iz)/(z), where /(z) was defined in p.39p . Since H is 
finite and because of p.40p . (A3) is now established. 

c) In order to conclude the proof of item 1., it remains to show that L is an (J^t)-martingale which is 
strongly orthogonal to M . This can be established similarly as in [53], Proposition 3.1, by making use of 
Fubini's theorem and Fubini's theorem for stochastic integrals (cf. (31], Theorem IV. 46) and (Al). (A2), (A3). 



Consequently, iHo, C , L ) provide a (possibly complex) FS decomposition of H. 
d) It remains to prove item 2., that is to say that the decomposition is real-valued. Let iHo,(^^ , L^) and 

^ 

(ifo,^ :L ) be two FS decomposition of H. Consequently, since H and (S't) are real- valued, we have 

= H-H =iHa -Ho)+ / iC" -lf )dSs + (L^ -l"), which implies that = /m(iJo) + Jo Irni^")dSs + 

Jq 

ImiLif). By Theorem [H] 1 . , the uniqueness of the real-valued Follmer-Schweizer decomposition yields that 
the processes (i/t),(^/^) and iLf) are real-valued. □ 



3.7 Representation of call and put options 



We used some integral representations of payoffs of the form p.27p . We refer to [T3], [32] and more re- 
cently [T7], for some characterizations of classes of functions which admit this kind of representation. In 
order to apply the results of this paper, we need explicit formulae for the complex measure 11 in some example 
of contingent claims. Let A' > be a strike. 



The European Call option H = iSx — -?'')+• For arbitrary < i? < 1, s > 0, we have 

(,s - K)+ 



1 

27ri 



-R+ioo fC^^^ 
R-ioo ziz — 1) 



(3.42) 
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The European Put option H = {K — St)+- For an arbitrary i? < 0, s > 0, we have 

{K - s)+ = — / s^- -dz . (3.43) 

4 The solution to the minimization problem 

FS decomposition will help to provide the solution to the global minimization problem. Let X be an additive 
process with cumulant generating function k. We denote 5*4 = soGxp(Xt),i € [0,T], sq > 0. Next theorem 
deals with the case where the payoff to hedge is given as a bilateral Laplace transform of the exponential of 
the additive process X. It is an extension of Theorem 3.3 of [53] to additive processes with no stationary 
increments. 

Theorem 4.1. Let H ~ /{St) where f is of the form VS. 27]) . We assume the validity of Assumptions\^\^ 
O The variance- optimal capital Vq and the variance- optimal hedging strategy Lp, solution of the minimization 
problem \2.2\ . are given by Vq — Hq and the implicit expression 



ft = + ^{Ht- -V^- f ipsdSs) , (4.1) 

Jo 



where the processes (Ht), (^t) and (At) are defined by 

j{z,t) := '^^ with pt{z,y) ^ Kt{z + y) - Kt{z) - Kt{y) , 

apt 

r]{z,dt) := Kdt{z) - -f{z,t)Kdt{l), At := ^^^^ 

dpt 

Ht := /e^'^''("^'^")5t"n(dz), -.^ f -f{z,t)e^*^ '^^'''^'^ S^I^Il{dz) . 
Jc Jc 

The optimal initial capital is unique. The optimal hedging strategy (pt(j-o) is unique up to some {¥{dLo) <Si dt)- 
null set. 

Remark 4.2. The mean variance trade-off process can be expressed as, see liS.ll]) . Kt — '^d"^^'* ^dtt(l)- 



Proof (of Theorem 14. ip . 

Since K is deterministic, the optimality follows from Theorem 13.301 and by items 1. and 2. of Theorem l2.7l 
We recall that a was given in p.lOp . Uniqueness follows from Theorem 12.61 2. □ 

When the underlying price is an exponential of additive process, we evaluate the so called variance of 
the hedging error of the contingent claim H i.e. the quantity E[(Vb + Grif) — H)"^], where Vb, and H 
were defined at Theorem |4JJ 



Theorem 4.3. Under the assumptions of Theorem ]^. 1\ the variance of the hedging error equals 

Jo-=( f I Mv,zmdy)Il{dz) 



where 



: otherwise, 
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with 

{—^ — ) '^^'^ ' 

P[y,z,t) := — . 4.2 

apt apt apt 

This expression of the error involving the function (3 (|4.2p , can be used to characterize the price models 
that are exponential of additive processes for which the market is complete, at least for vanilla option payoffs. 
For instance, by evaluating /?, we can verify, in Remarks 15.101 and 15.111 below, the complete market model 
property in the Poisson and the Gaussian case. 

Proof (of Theorem 14. 3[) . Since Xq = 0, is the trivial cr-field, therefore Lq = 0, because it is mean-zero 
and deterministic. 

The quadratic error can be calculated using Theorem 12.71 3. It gives 

fT 



E 



/ exp{~{KT-Ks)}d{L'') 
Jo 



(4.3) 



where is the remainder martingale in the FS decomposition of H. We proceed now to the evaluation 
of (L^)- Similarly to the proof of Theorem 3.2 of [M], using p.23p . the bilinearity and the stability w.r.t. 
complex conjugate of the covariation together with p.4ip , it is possible to show that 



{L^,L")^ =J J {L{y),L{z))^n{dy)Ilidz). (4.4) 
It remains to evaluate {L{y), L{z)) for y,z & supp(n). We know by Proposition 13 . 1 71 that for all z € -jj 

(M(y),M(z)),= f' Slli'pduiy.z) . 



Using the same terminology as in Proposition 13. 25i similarly to (|3.24[) we have 
{L{y), L{z)), = /* eijM^,du)+,iy4u))sy+^ [^^^(y^ ^) _ ^(^^ 



t 



We come back to (|4.3p . Recalling that a{y, z, t) ~ {ri{z, T) — ri{z, t)) — {r]{y, T) — rj{y, t)) — [Kt — Kt), where 
K is the MVT process, we have 

-(^T-K.)d {L{y), Liz)), = r e"('''^'*)5j'+^/3(y, z, t)dpt. 

Jo 

Since E[iS'fjl'^] = SQ~^^e'^*^y^^\ an application of Fubini's theorem yields 

E e-(^--^')d (L(y), L{z))^ - £ ^.(y,z,t)+.,iy+.) ^^y^ ^^^^^^ (4 5) 

which equals Jo(jJ, z). (|4.4p . (|4.5p and again Fubini's theorem imply 

^HKt-k,)^ (i^, L"), ^ j^j^ My, z)n{dy)U{dz). 
This concludes the proof of Theorem 14.31 □ 
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5 Examples 



5.1 Exponential of a Wiener integral driven by a Levy process 

Let A be a square integrable Levy process and let {t,z) i-> n^{z) be the cumulative generating function of 
A with domain Da in the sense of Definition 12.11 (i, z) H' n^{z) is continuous because of Proposition 13.251 
We observe that 

4{z) = tK^{z) , (5.6) 

where : A ^ C is a continuous function such that k.^{z) = Hi[z). Let Z : [0, T] — > M be a bounded 
Borel function. We will consider in this subsection the additive process Xt = Jq IsdAs- Let us define the set 
L'a(0 C M such that 

Da{1) ^ {x eR\ixJx e Da} + iR , where ^ := inf /, I:=sup/. 

Lemma 5.1. The cumulant generating function of X is such that for all z e Da{1), we have 

KXt{z) = / ti^{zls)ds. 
Jo 

In particular D\(l) C D, where D is the domain defined according to Definition \2.1\ 

Proof. If I is continuous, the result follows from the observation that jj" IgdAs is the limit in probability of 
Z]^=o ^tj (^tj+i ~ ^tj) where ^ to < ti < ... < tp = T is a subdivision of [0,r] whose mesh converges to 
zero. If I is only Borel bounded the result can be established through approximation by convolution. □ 

We formulate the following hypothesis which will be in force for the whole subsection. 

Assumption 4. 7. k^(2) - 2k^(1) 7^ 0. 

2. l_>0 and 21 e Da. 

Remark 5.2. Lemma \3.9\ applied to X being the Levy process A implies that, for every 7 > 0, such that 
27 G Da, we have 

K^i2-f) ~ 2k^(7) > . (5.7) 

Remark 5.3. 1. By item 2. of Assumption^ 2 e Da{1) and so does 1 because Da{1) is convex. By 
Lemma \5.1[ 1 and 2 belong to D. 

2. pt = J^{K^i2ls)~2fi^ih))ds; 

3. X is a semimartingale additive process since t 1—^ Kt(2) has bounded variation, see Proposition \3.4\ 
Proposition 5.4. Assumptions{I\ and\^ are verified. Moreover Da{1) C V. 

Proof. 1. By item 1. of Remark 15.31 2 e D and so the second item of Assumption [1] is verified. Using 
Lemma [3.91 item 1. of Assumption [T] is verified if we show that t i-> pt{l) — Kt(2) — 2Kt(l) is strictly 
mcreasmg. Now Kt(2)-2Kt(l) {k^{21s) - 2k^{Is)) ds. Inequality and item 2. of Assumption 
pimply that Vs G [0,r]. k^{21s) — 2k^{Is) > 0, and consequently t i-> /Ot(l) is strictly increasing. 
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2. For z G D\{1), by Lemma [5TT] and Remark 15.31 2. we have 



dKt{z) 




dpt 





Previous supremum and infimum exist since x i— > k^{zx) is continuous and it attains a maximum and 
a minimum on a compact interval. So, -Da(Z) C 2? and Assumption [2] is verified because of point 1. in 
Remark 15.31 

□ 

Remark 5.5. Suppose for a moment that 

2IoC{x\lx,lxG Djy}. (5.9) 

1. That implies then 2Iq C Da(0- Point 3. of Assumption\^ follows by Lemma \5.1\ Item 2. of the same 
Assumption is also verified. In fact, since 2Iq C D\{1) and 2 e D\{1) and because of the fact that 
D\{1) is convex, we have 

suppn u (suppn + 1) c ^ + ^ c D^{i). 

The conclusion follows by Proposition \5.4\ which says Da(0 *^ ^• 

2. From the proof of Proposition \5.4\ it follows that 

dpt ~ K^{2lt) - 2n^{lt) ■ 
Admitting point 1. of Assumption\^ then I is compact. Taking into account (j5.8p . the fact that 
1 G D\{1), so I C D\{1), and that is continuous, point 4- of Assumption[3[is verified. 

We consider again the same class of options as in previous subsections. To conclude the verification of 
Assumption 13] it remains to show the following. 

• Iq is compact. This point will be trivially fulfilled in the specific cases. 

• (EH). 

Example 5.6. We keep in mind the call and put representations provided in Section \3. 71 

1. H ^ {St — K)+- In this case 2Iq = {2R, 2} and (|5.9p is verified, since R g]0, 1[. 

2. H = {K - St)+- Again, here R<0, 2Io = {2R}. 

Again, we only have to require that D\ contains some negative values, which is the case for the three 
examples introduced in Remark \5.8[ Selecting R in a proper way, (j5.9p is fulfilled. 

Corollary 5.7. We consider a process X of the form Xt — /J IsdAg under Assumption^ The FS decom- 
position of an option H of the type (j3.27p and the related solution to the minimization problem are provided 
by Theorem \3.3(A Provosition \3.25\ and Theorem \4.1\ together with the expressions given below. 
For z G suppll, t G [0, T] we have 

i^^'jls) , _ n^{{z + l)h)-K''{zls)-K''{ls) 

k^{21s)-2k^{Is)' k^{21s)-2k^{Is) 

V{z,s) = n\zh) - ^f,^2lt-2nHls) ^^^^^^ ^ ^^^'^ ~ """^^"^'^ ' """^^^'^^ • 
Again, for convenience, if z ^ suppII then we define 7(2, •) = 77(2, •) = 0. 
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5.2 Considerations about the Levy case 

li I = 1 then X coincides with the Levy process A and Assumption 2] is equivalent to Hubalek et aha Con- 
dition introduced in [M] i.e. 1. 2e D; 2. ^^^(2) - 2k^(1) ^ . 

In that case we have D = D\ ~ D\(l). Therefore V ~ D because ^^(2) = . . for 

dpt (2) — 2k'^(1) 

any t e [0,r],z e D. 

We recall some cumulant and log-characteristic functions of some typical Levy processes. 

Remark 5.8. 1. Poisson Case: If X is a Poisson process with intensity X, we have that k^{z) = A(e^ — 1). 
Moreover, in this case the set D\ — C. 



2. NIG Case: This process was introduced by Barndorff- Niels en in J^. Then X is a Levy process with 
Xi ^ NIG{a,l3,S,fi), with a > > 0, (5 > and p £ R. We have k^{z) = fiz + S{jo — jz) and 
-/z = ^/d^~[^Tzf, Da = [-a- /3,a- f3]+iR. 



3. Variance Gamma case: Let a, /3 > 0, (5 7^ 0. If X is a Variance Gamma process with Xi ^ VG{a, /3, 5, fi) 
with K^{z) ^ fiz + 6Log I — — ^ ) , where Log is again the principal value complex logarithm de- 



fined in Section 2. The expression of k (z) can be found in \24\ \27j or also flS^ . table IV. 4-5 



m 



the particular case /i = 0. In particular an easy calculation shows that we need z € <C such that 
Re{z) e] - f3 - y//3'^ + 2a, -~f3 + y/ f3'^ + 2a[ so that k^{z) is well-defined so that 



Da=]-I3- ^/l3^ + 2a, -~(3 + ^/S^ + 2a[+iR. 

Remark 5.9. We come back to the examples introduced in Remark \5.8\ In all the three cases, Hubalek et 
alia Condition is verified if 2 £ D . This happens in the following situations: 

1. always in the Poisson case; 

2. if A = X is a NIG process and if 2 < a — (3 ; 



3. if K^X IS aVG process and i/2 < + ^13"^ + 2a . 
Theorem 14.11 allows to re-obtain the results stated in |24| . 
Remark 5.10. If X is a Poisson process with parameter A > then the quadratic error is zero. In fact, 

K^{z) = \{exp{z) - 1)) , pt{y, z) = \t{exp{y) - l)(exp(2:) - 1) 

K^{z + 1)-K^{z)-K^{l) _exp(z)-l 
^^"^'^^ ~ kA(2)-2kA(i) e-1 

imply that f3{y, z,t) ~ for every y, z £ C,t £ [0, T] . 

Therefore Jo{y, z, t) = 0. In particular all the options of type ^3.27]) are perfectly hedgeable. 
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5.3 About some singular non- stationary models 

Here, we consider some singular models, in the sense that the cumulant generating function of the log-price 
process is not absolutely continuous with respect to (a.c. w.r.t.) Lebesgue measure. More precisely, let (Wt) 
be a standard Brownian motion. A classical approach to model the volatility clustering effect consists in 
introducing the notion of trading time (as opposed to the real time) which accelerates or slows down the 
price process depending on the activity on the market. This virtual time is represented by a change of time 
iTt)t>a and the log-price is then constructed by subordination i.e. Xt = W^^ty Now, if the change of time r 
is singular, then it can be proved that the log-price process X is also singular. 

This typically happens when the change of time r, is obtained as the cumulative distribution function of a 
deterministic positive niultifractal measure dT{t) = dip(t), singular w.r.t. Lebesgue measure. Multifractal 
measures were introduced in the physical sciences to model turbulent flows [28]. More recently, in [10], the 
authors used this construction precisely for modeling financial volatility. But their model, the Multifractal 
Model of Asset Returns (MMAR), relies on a random (and not deterministic) multifractal measure and is 
hence beyond the framework of this paper. 

Below, we consider two examples of singular non-stationary log-price models based on such (deterministic 
or random) singular changes of time. 

1. Deterministic change of time (log-Gaussian continuous process): Let us consider the log-price process 
X such that Xt = W.^(^t) ■ where ip : K+ M+ is a strictly increasing function, including the pathological 
case where ■^jJ^ = a.e. For z G D = C, we have £[6"^-^'] = E[e^'^'''(-)] = e'^*'^"^ = e^"^*^*), so that 
^((z) = ^il){t) , pt = ipit). Notice that dKt{z) is not necessarily a. c. w.r.t. Lebesgue measure and 
that this is verified as soon as d'ipit) <€i dt. Assumption [Tj 1. is verified since ip is strictly increasing; 
Assumption [1] 2., Assumption [2] and Assumption [3] are verified since D = V = C and ^^^^^ — ^ 
continuous. Consequently all the conditions to apply Theorem 14. II are satisfied and 

/ \ / \ '0(^) / 2 \ 1 

7(z,t)=z, Tj{z,t} = —^{z - z) and At = -. 

Remark 5.11. Calculating j3{y,z,t) in (|4.2p . we find P = 0. Therefore here also the quadratic error 
is zero. This confirms the fact that the market is complete, at least for the considered class of options. 

2. Random change of time: Let {9t)t>o denote an increasing Levy process such that 9i follows an Inverse 
Gaussian distribution with parameters 5 > and 7 > 0. Now, let us consider Y the process such that 
Yt ^ lit + I39{t) + We(t). for all t G [0, T], with ;3 , ^ e R. Then one can prove that y is a NIG Levy 
process with Yi ^ NIG{a = \/7^ -f- /3^, /3, 5, fi). Finally, let us consider the log- price process X such 
that Xt = Wrt, where Tt = 9ii>{t) a-nd tp is the cumulative distribution of a deterministic multifractal 
measure on [0,r]. Hence, the cumulant generating function of Xt is singular w.r.t. Lebesgue measure 
and is given by Kt{z) = [iiz + S{'^o ^ 7z)]V'(0 with 7z = -^/a^ — (/? + 2)^, for all z £ D := Dxt = 
[-a- P,a- P] +iR. 

6 Application to Electricity 

6.1 Hedging electricity derivatives with forward contracts 

Because of non-storability of electricity, no dynamic hedging strategy can be performed on the spot market. 
Hedging instruments for electricity derivatives are then futures or forward contracts. For simplicity, we will 
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assume that interest rates are deterministic and zero so that futures prices are equivalent to forward prices. 
The value of a forward contract offering the fixed price F^^'' at time for delivery of IMWh at time Td is by 
definition of the forward price, Sq''^'' = 0. Indeed, there is no cost to enter at time the forward contract 
with the current market forward price Fq''. Then, the value of the same forward contract 5°'"^'* at time 
t e [0,Td\ is deduced by an argument of Absence of (static) Arbitrage as 3°''^'' = e"'''^^''"*) (F/''' - F^''). 
Hence, the dynamics of the hedging instrument {S^''^'')o<t<Ta is directly related (for deterministic interest 
rates) to the dynamics of forward prices {F^'^)(xt<Td- Consequently, in the sequel, when considering hedging 
on electricity markets, we will always suppose that the underlying is a forward contract (5'('^'^'*)o<t<Td and 
we will focus on the dynamics of forward prices. 

6.2 Electricity price models for pricing and hedging application 

Observing market data, one can notice two main stylized features of electricity forward prices: 

• Volatility term structure of forward prices: the volatility increases when the time to maturity decreases. 
Indeed, when the delivery date approaches, the flow of relevant information affecting the balance 
between electricity supply and demand increases and causes great variations in the forward prices. 
This maturity effect is usually referred to as the Samuelson hypothesis, it was first studied in |34| and 
can be observed on Figure [TJ in the case of electricity futures prices. 

• Non-Gaussianity of log-returns: log-returns can be considered as Gaussian for long-term contracts but 
begin to show heavy tails for short-term contracts. 




1 2 3 4 5 6 

Time to maturity (in month) 



Figure 1: Volatility of electricity Month-ahead futures prices w.r.t. to the time to maturity estimated on the French 
Power market in 2007. 

Hence, a challenge is to be able to describe with a single model, both the non-Gaussianity on the short term 
and the volatility term structure of the forward curve. One reasonable attempt to do so is to consider the 
exponential Levy factor model, proposed in [3] or |I2| . The forward price given at time t for delivery at time 
Td > t, denoted F^ is then modeled by a p-factors model, such that 

p 

F^" = F^" cxpimj" + Xt'^" ) , for aU t e [0, Td] , where (6.10) 
fc=i 
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• {mJ'^)o<:t<Ta is a real deterministic trend; 

• for any fc = 1, • • -p, {Xt'^'')o<t<T^ is such that X^'^" = /J crfce-^'=(^''-")dA^\ where A = (A\ • • • , A^) 
is a Levy process on K'', with E[A5^] = and ya7'[A5^] = 1: 

• CTfc > , Afc > , are called respectively the volatilities and the mean-reverting rates. 

Hence, forward prices are given as exponentials of additive processes with non- stationary increments. In 
practice, we consider the case of a one or a two factors model (p = 1 or 2), where the first factor is a 
non-Gaussian additive process and the second factor is a Brownian motion with ci ^ a2. Notice that 
this kind of model was originally developed and studied in details for interest rates in |32| , as an extension of 
the Heath- Jarrow-Morton model where the Brownian motion has been replaced by a general Levy process. 
Of course, this modeling procedure (|6.10|) . implies incompleteness of the market. Hence, if we aim at pricing 
and hedging a European call on a forward with maturity T < T^, it won't be possible, in general, to hedge 
perfectly the payoff {Fr^'' — K)+ with a hedging portfolio of forward contracts. Then, a natural approach 
could consist in looking for the variance optimal initial capital and hedging portfolio. In this framework, the 
results of Section |3] generalizing the results of Hubalek & al in [24] to the case of non stationary additive 
process can be useful. 

6.3 The non Gaussian two factors model 

To simplify let us forget the superscript Td denoting the delivery period (since we will consider a fixed 
delivery period) . We suppose that the forward price F follows the two factors model 

Ft = Fo exp(TOt + Xl + Xf) , for all t e [0, T<j] , where (6.11) 

• m is a real deterministic trend starting at 0. It is supposed to be absolutely continuous w.r.t. Lebesgue; 

• X^ = crse-^(^'i-")dA„, where A is a Levy process on M with A following a Normal Inverse Gaussian 
(NIG) distribution or a Variance Gamma (VG) distribution. Moreover, we will assume that E[Ai] = 
and yar[Ai] = 1; 

• X^ = aiW where W is & standard Brownian motion on R; 

• A and W are independent; 

• ffs and ai standing respectively for the short-term volatility and long-term volatility. 

6.4 Verification of the assumptions 

The result below helps to extend Theorem 14.11 to the case where X is a finite sum of independent semi- 
martingale additive processes, each one verifying Assumptions [1] [2] and [3] for a given payoff H = f{soe^^). 

Lemma 6.1. Let X^,X^ be two independent semimartingale additive processes with cumulant generating 
functions and related domains D'^V^ ,i ~ 1,2 characterized in Remark \2.1U and (|3.12p . Let / : C — > C of 
the form IKFA) . 

For X = X^ + X^ with related domains D,T> and cumulant generating function k, we have the following. 
1. D = D^n D'\ 
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2. 2?i n 2?2 c V. 

3. If X^,X^ verify Assumptionsl^\^ and\^ then X has the same property. 

Proof. Since X^,X^ are independent and taking into account Remark 12.101 we obtain 1. and Kt{z) = 
+ K^{z), Vz e D. We denote by = 1,2, the reference variance measures defined in Remark 13.121 
Clearly p = + and dp^ <C dp with 
If z G n I?^. we can write 



dp* I 



< 1. 



dKt{z) 



dpt 



dpt 



< 



dnjiz) dpi 



dpi dpt 



dpt 



^ di4{z)dpl 



dpi dpt 



dpt 



duKz) 



< 



dpi 
dnliz) 



'dpi 



dpi 



dpi 
dpt 

2 

dpi 



dnKz) 



dpi 



dpi 

dpt 



dp 



di4{z) 



dpi 



dpi 



This concludes the proof of I?^ n I?^ C I? and therefore of the of Point 2. 
Finally Point 3. follows then by inspection. 



□ 



With the two factors model, the forward price F is then given as the exponential of an additive process, 
X, such that for all t e [0, Td], 

rt 

e'^^^'-'^UAu + aiWt . (6.12) 







Xt=mt+ Xl + Xl ^mt+<Ts 

For this model, we formulate the following assumption. 
Assumption 5. 1. 2a s E Da. 

2. If (7i = 0, we require A not to have deterministic increments. 

3. f IS of the type fulfilling ^ 



Proposition 6.2. 1. The cumulant generating function of X defined by \6.12\) . k : [Q,Td] x D 
such that for all z G Dj^[as) and for all t G [0, T^], 



z^aft 



K,(z) = zmt + ^^+ / «:^(z(T,e-^(^^-"))dM . 



(6.13) 



In particular for fixed z £ Djy{as), t H- Kt(z) is absolutely continuous w.r.t. Lebesgue measure. 

2. Under Assumptions^ Assumptions^^ \^and\^ are fulfilled. 

Proof. We set X"^ = m + X'^. We observe that = = C, (z) — cxp{zmt + z^a'f^). We recall that 
A and W are independent so that X^ and X^ are independent. For clarity, we only write the proof under 
the hypothesis that A has no deterministic increments, the general case could be easily adapted. X^ is a 
process of the type studied at Section [5TT1 it verifies Assumption [4] and D\{1) contains Da{o's). 

According to Proposition 15. 4i Remark 15.51 and (|5.9p it follows that Assumptions [1] [2] and [3] are verified 
for X-^. Both statements 1. and 2. are now a consequence of Lemma WA\ □ 



The solution to the mean- variance problem is provided by Theorem 14.1 
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Theorem 6.3. We suppose Assumptions^ The variance- optimal capital Vo and the variance-optimal hedg- 
ing strategy ip, solution of the minimization problem (|2.2I) , are given by Theorem \4-l\ and Theorem \ cl.3(A 
Proposition \3. 25\ together with the expressions given below: 



If. = 

-f{z,t): = 

T]{z,t): = 



zaf + k\{z + l)lt) - K^izk) - K^il)t 



af + K^{2lt) - 2t,^{lt) 

-A/ 



zmt + + ^'^{zlt) - -i{z,t){nn + Y + 



dt 



Remark 6.4. Previous formulae are practically exploitable numerically. The last condition to be checked is 

2g, G Dt^. (6.14) 
1. Ai is a Normal Inverse Gaussian random variable; if as < then (j6.14p is verified. 



2. Ai is a Variance Gamma random variable then (|6.14p is verified; if for instance as < ^^^^ ~~ 



7 Simulations 

We are interested in comparing, in simulations, the Variance Optimal (VO) strategy to the Black-Scholes 
(BS) strategy when hedging a European call, with payoff {St — -f^)+, on an underlying stock with log-prices 
Xt = log(S't) that have independent but non Gaussian increments. More precisely, we assume that the 
underlying is an electricity forward contract St — S^''^'^ = e~^^'^'^~*'\Fj''^ — Fq'') with delivery date Td equal 
to the maturity of the call Td = T. 

First, we consider the case where the log-price process X is an exponential of a Levy process, continuing the 
analysis of [24] , then we consider the non stationary case. We make use of different simulated data according 
to the underlying model, stationary in one case, non stationary in the second one. 

Our simulations investigate two features which were not considered in |24) (even in the stationary case): 
first the robustness of the BS hedging strategy w.r.t. the underlying price model, second the sensitivity of 
the continuous VO strategy w.r.t. to the discreteness of the trading dates. 

The VO strategy knows the real incomplete price model (with the real values of parameters) whereas the 
BS strategy assumes (wrongly) a log- normal price model (with the real values of mean and variance). Of 
course, the VO strategy is by definition optimal, w.r.t. the quadratic norm. However, both strategies (VO 
and BS) are implemented in discrete time, hence our goal is precisely to analyze the hedging error outside of 
the theoretical framework of a continuously rebalanced portfolio. Moreover, we are interested in interpreting 
quantitatively the differences between both strategies w.r.t. to some characteristics such as the underlying 
log-returns distribution or the number of trading dates. 

The time unit is the year and the interest rate is zero in all our simulations. The initial value of the underlying 
is sq = 100 Euros. The maturity of the option is T = 0.25 i.e. three months from now. 



7.1 Exponential Levy 

In this subsection, we simulate the log-price process X as a NIG Levy process with Xi ^ NIG{a, /3,S, fi). 
Five different sets of parameters for the NIG distribution have been considered, going from the case of 
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almost Gaussian returns corresponding to standard equities, to the case of highly non Gaussian returns. 
The standard set of parameters is estimated on the Month-ahead base forward prices of the French Power 
market in 2007: 

a = 38.46 , ^ = -3.85 , S = 6.40 , = 0.64 . (7.15) 

Those parameters imply a zero mean, a standard deviation of 41%, a skewness (measuring the asymmetry) 
of —0.02 and an excess kurtosis (measuring the fatness of the tails) of 0.01. The other sets of parameters are 
obtained by multiplying the parameter a by a coefHcient C, (/3, S, fi) being such that the first three moments 
are unchanged. Note that when C grows to infinity the tails of the NIG distribution get closer to the tails of 
the Gaussian distribution. For instance, Table [T] shows how the excess kurtosis (which is zero for a Gaussian 
distribution) is modified with the five values of C chosen in our simulations. 



Coefficient 


C = 0.08 


C = 0.14 


C = 0.2 


C = 1 


C = 2 


a 


3.08 


5.38 


7.69 


38.46 


76.92 


Excess kurtosis 


1.87 


0.61 


0.30 


0.01 


4.10^3 



Table 1: Excess kurtosis of Xi for different values of a, [13,5, ji) insuring the same three first moments. 
7.1.1 Strike impact on the initial capital and the hedging ratio 

Figure [2] shows the initial capital (on the left graph) and the initial hedge ratio (on the right graph) produced 
by the VO and the BS strategies as functions of the strike, for three different sets of parameters C = 
0.08 , C = 1 , C = 2. We consider = 12 trading dates, which corresponds to operational practices on 
electricity markets, for an option expiring in three months. One can observe that BS results are very similar 
to VO results for C > 1 i.e. for almost Gaussian returns. However, for small values of C, for C ~ 0.08, 
corresponding to highly non Gaussian returns, BS approach under-estimates out-of-the-money options and 
over-estimates at-the-money options (for K = 99 Euros the BS initial capital is equal to 8.65 Euros i.e. 122% 
of the VO initial capital, while for K = 150, it vanishes to 23 Cents i.e. only 57% of the VO initial capital). 




strike price Strike price 



Figure 2: Initial capital (on the left) and hedge ratio (on the right) w.r.t. the strike, for C = 0.08 , C = 1 , C = 2. 
7.1.2 Hedging error and number of trading dates 

Figure [3] considers the hedging error (the difference between the terminal value of the hedging portfolio and 
the payoff) w.r.t. the number of trading dates, for a strike K =^ 99 Euros (at the money) and for five different 
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sets of parameters C given on Table [TJ The bias (on the left graph) and standard deviation (on the right 
graph) of the hedging error have been estimated by Monte Carlo method on 5000 runs. Note that we could 
have used the formula stated in Theorem 14.31 to compute the variance of the error, but this would have 
given us the limiting error which does not take into account the additional error due to the finite number of 
trading dates. 

In terms of standard deviation, the VO strategy seems to outperform noticeably the BS strategy, for small 
values of C ( for C = 0.08 the VO strategy allows to reduce 10% of the standard deviation of the error). 
As expected, one can observe that the VO error converges to the BS error when C increases. This is due 
to the convergence of NIG log-returns to Gaussian log-returns when C increases (recall that the simulated 
log-returns are almost symmetric). On Figure [3l the hedging error (both for BS and VO) decreases with the 
number of trading dates and seems to converge to a limiting error. Here, it is interesting to distinguish two 
sources of incompleteness, the rebalancing error due to the finite number of trading dates and the intrinsic 
error due to the price model incompleteness. For instance, one can observe that for small values of C < 0.2, 
even for small numbers of trading dates, the intrinsic error seems to be predominant so that it seems useless 
to increase the number of trading dates over > 12 trading dates. Moreover, surprisingly one can observe 
that for a small number of trading dates A < 12 and for large values of C > 1, BS seems to outperform 
the VO strategy, in terms of standard deviation. This can be interpreted as a consequence of the central 
limit theorem. Indeed, when the time between two trading dates increases the corresponding increments of 
the Levy process converge to a Gaussian variable. Similarly to the observation of [T^], section 5., in term 
of hedging errors, BS strategy seems to be quite close to VO strategy. The same kind of conclusions were 
obtained in the discrete time setting by [1]. 

In term of bias, the over-estimation of at-the-money options (observed for C ~ 0.08, on Figures [5]) seems 
to induce a positive bias for the BS error (see Figure [3]), whereas the bias of the VO error is negligible (as 
expected from the theory). 




■ 10 20 30 40 50 60 70 10 20 30 40 50 60 70 

Number of trading dates: N Number of trading dates: N 

Figure 3: Hedging error w.r.t. the number of trading dates for different values of C and for K — 99 Euros (bias, on 
the left and standard deviation, on the right). 



7.2 Exponential of additive processes 

In this subsection, we simulate the log-price process X as an additive process such that 

Xf = asC"^^^'^^ dKu where A is a Levy process with Ai ^ NIG{a, /3^S, fi) 
Jo 
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The standard set of parameters (C = 1) for the distribution of Ai is estimated on the same data as in the 
previous section (Month-ahead base forward prices of the French Power market in 2007): 

a = 15.81 , P = -1.581 , 6 = 15.57 , /i = 1.56 . 

Those parameters correspond to a standard and centered NIG distribution with a skewness of —0.019. The 
estimated annual short-term volatihty and mean-reverting rate are cr^ = 57.47% and A = 3. The other sets 
of parameters considered in simulations are obtained by multiplying parameter a by a coefficient C, S, fi 
being such that the first three moments are unchanged). 

The results are comparable to those obtained in the case of the Levy process, on Figure 01 




°0 5 10 15 20 25 30 35 40 45 50 \ 5 10 15 20 25 30 35 40 45 50 

Number of trading rlates: N Number of trariing riafes: N 

Figure 4: Hedging error w.r.t. the number of trading dates for C = 0.08 and C = 1, for K = 'd'd Euros (bias, on the 
left and standard deviation, on the right). 
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